Crossover and scaling in a two-dimensional field-tuned superconductor 
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Using an analysis similar to that of Imry and Wortis, it is shown that the apparent first order 
superconductor to metal transition, which has been claimed to exist at low values of the magnetic 
field in a two-dimensional field-tuned system at zero temperature, can be consistently interpreted 
as a sharp crossover from a strong superconductor to an inhomogeneous state, which is a weak 
superconductor. The true zero-temperature superconductor to insulator transition within the inho- 
mogenous state is conjectured to be that of randomly diluted XY model. An explaination of the 
observed finite temperature approximate scaling of resistivity close to the critical point is speculated 
within this model. 
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Two-dimensional (2D) zero-temperature magnetic 
field(i/ )-tuned superconductor-insulator transition (SIT) 
has recently attracted a lot of attention [1-7]. It is com- 
mon to think of this transition as a prototypical contin- 
uous quantum phase transition (QPT) at a critical mag- 
netic field H = H c with a diverging correlation length 
£ ~ \H C — H\~ v , and associated scaling [8]. Here, v is the 
correlation length exponent. The state on the H < H c 
side of the transition is a superconductor where the vor- 
tices are localized in a so-called vortex glass state, and 
the Cooper pairs are delocalized. The state on the other 
side of the transition is an insulator, where the vortices 
form a superfluid (they follow Bose statistics as well in 
a dual representation [9,10]), and the Cooper pairs are 
localized. Exactly at the transition, both the vortices 
and the Cooper pairs are delocalized, and the system is 
supposed to show metallic behavior with a universal re- 
sistivity [8,11,12]. 

However, recent experiments on amorphous Mo-Ge 
films do not show scaling in the limit of very low tem- 
peratures T < lOOmK [4], do not show universal re- 
sistivity at the putative transition point, and, in fact, 
show signatures of zero temperature saturation of Ohmic 
resistance in a wide range of H on both sides of the 
putative H c [1,4]. The system thus appears to have a 
metallic ground state in almost all cases except below 
a very low value of the magnetic field, H sm ('sm' for 
superconductor-metal), where true superconductivity is 
apparently recovered [6]. Hysteresis is observed in the 
vicinity of H sm . Phenomenologically, there seems to be a 
zero-temperature first order field-tuned superconductor 
to metal transition at H ~ H sm . The conventional pic- 
ture of the field-tuned SIT, then, needs to be seriously re- 
visited. There is also experimental evidence of a metallic- 
like phase in thickness or disorder-tuned superconductor- 
insulator systems [13,14]. The possibility of a metallic 
phase intervening between the insulating and the super- 
conducting phases has been a subject of many recent the- 
oretical papers [3,15-17]. In the following, we will discuss 
only the field-tuned systems; whether similar considera- 
tions can be applied to other systems as well is left for 



future study. 

We will first focus on the low field superconductor to 
metal transition in the field-tuned system and argue that 
it is a remnant effect of the effects of quenched disorder 
on the first order vortex lattice melting transition of the 
pure system at a higher value of H — H cr . Effects of 
quenched disorder on first order transitions have a long 
history in the literature. Imry and Ma [18] showed that 
for disorder coupling as random fields conjugate to the 
order parameter, first order transitions are completely 
eliminated below a critical dimension which depends on 
the symmetry of the order parameter. Imry and Wortis 
[19] then argued that for disorder coupling as random 
bond (or random T c in a classical system), the disconti- 
nuities associated with a first order transition are elim- 
inated. They also proposed a mechanism of how that 
happens and indicated the possibility that for this kind 
of disorder-coupling the first order transition may re- 
duce to a continuous one. Later, Aizcnman and Wehr 
[20] proved mathematical theorems conclusively proving 
these heuristic assertions. In this paper, we will closely 
follow the argument of [19], but make the most use of a 
feature, often overlooked, that the elimination of discon- 
tinuities of a first order transition due to random-bond 
coupling of the disorder can start as a sharp crossover 
at D = 2. Beyond the sharp crossover and within the 
resulting inhomogenous state, the true thermodynamic 
phase transition can be percolative, which is a continu- 
ous transition and thus consistent with the established 
theorems. 

We will work in terms of field-induced vortices instead 
of Cooper pairs. The putative SIT is then the melting 
transition of the vortex solid to a vortex liquid, which, in 
fact, is a vortex superfluid [8]. In the limit of a clean sys- 
tem free of any impurities, this is the first order Abrikosov 
lattice melting transition at n = n CI . For an average 
impurity concentration p, let's assume that n cr is renor- 
malized to n cr (p). At this level, n cv (p) is an assumed 
first order transition point between a vortex solid phase 
and a vortex liquid phase. The vortex solid phase may 
be a vortex glass (see, however, [21]), or a phase where 
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the vortices are trapped individually by the impurity po- 
tentials, or even some (as yet) unknown state where the 
vortices are localized (note that in the presence of ran- 
dom pinning by impurities, which couple to a crystalline 
solid as random fields, the translational long range or- 
der of the crystal is destroyed below D = 4 [18]). Our 
second assumption is that, in the vortex solid phase, 
local fluctuations in the impurity concentration modu- 
late the local n cr . This is sensible (although we cannot 
prove it here), since the impurities, along with the re- 
pulsive interaction among the vortices, are responsible 
for vortex-localization. This can also be seen as follow- 
ing from the dependence of n cr (p) on p. Local p may 
fluctuate from the average p, and locally modulate the 
critical concentration. Later, we will argue, following 
[19], that n cr (p) will be rounded in a real system by do- 
main formation on both sides of the transition, and the 
hypothetical vortex solid will loose its meaning. This is 
the reason why a precise characterization of the vortex 
solid is irrelevant for the present purpose. This line of 
argument will, however, bring out the importance of a 
sharp crossover-concentration, n sm , which we will iden- 
tify with H sm , beyond which the domains proliferate and 
the system shows appreciable low-temperature resistance 
[3] . The true macroscopic superconductor-insulator tran- 
sition at H c can then only be defined in terms of perco- 
lation, by that of randomly diluted Josephson Junctions 
between the domains of the vortex solid phases. Since the 
zero-temperature SIT in this framework is essentially ge- 
ometric, and not a prototypical QPT, scaling at finite 
temperature close to the critical point is not immediate. 
We will speculate on the consequences of such a transi- 
tion to finite temperature scaling of Ohmic resistivity. 

In order to clarify the said remnant effect, let us flesh 
out the arguments for the effects of disorder on a first 
order transition. In our formulation in terms of a hypo- 
thetical vortex solid (and not a crystalline lattice), the 
analysis of [19] directly applies. It will pay to briefly go 
through the argument, however, since we will use the fi- 
nal result to explain the experiments. It will also help 
develop the alternative framework of SIT presented here 
(see also [3]). Since the correlation length £ cr remains 
finite at the transition, only regions of finite area ~ 
are correlated close to the transition. The typical den- 
sity fluctuation of the impurities in a correlated region of 
linear dimension £ cr is [19] 



Ap~[p(l-p)]'G 1 - 



(1) 



Density fluctuations of other magnitudes are possible, 
but their probabilities are exponentially suppressed [22]. 
So a correlated region has impurity concentration p + 
Ap > p > p — Ap with maximum probability. Now sup- 
pose we are approaching the transition point by increas- 
ing the perpendicular magnetic field H, thus increasing n. 
As long as n < n cr (p) the system is nominally in the solid 



phase. But a typical correlated region may already be in 
the wrong phase due to fluctuations in the impurity con- 
centration. The wrong phase is a vortex liquid, which is 
not yet a supcrfluid since the size is finite. Thus, if there 
were no interfacial surface energy between the solid and 
the liquid phases, then the first order transition would 
have had a smearing width 



dn cr (p) 
r = I dp \ A P> 



(2) 



and the domain-sizes would have been of the order of £ cr . 

But if we take into account the interfacial surface en- 
ergies, as we must, then the domain-sizes can be much 
bigger than £ cr . This is possible, because the system has 
to balance the loss in the surface energy by gaining suf- 
ficient volume energy. By creating domains comprising 
several correlated regions with linear dimensions £ cr , it 
can maximize the gain in volume energy. One can show 
that [19], as one approaches the critical concentration 
n CI (p) by varying H, clusters of a general size I will be 
first unstable at a concentration n cr (p) — An whenever 



An < 5(0, 

where g(l) is a function of the domain size I, 



(3) 



g(l) = \ P( 1-P)]^\ d ^\ - gtTn :ff r2 . (4) 



dp 



Here a is the interfacial surface tension, C is a geo- 
metrical factor of order one, and f(p) is a function of p, 

f(p) = [" a(/ 3n /3) ]p,"c,(p)' where ft and f s are the free 
energy densities of the vortex liquid and solid phases re- 
spectively. The formula in Eq. 4 is for D = 2. The first 
term is a measure of the volume energy gain by creating 
the liquid phase in a domain of size I, and the second term 
comes from the surafce energy lost by doing it. Equation 
3 is the first order analog of the Harris criterion [23]. 
Both the rounding width An and the domain-size /* are 
now determined by the single equation, 



An = g(l*), 



(5) 



where /* maximizes g(l). This does not mean, however, 
that for n < n CI (p) — g(l*) the system is free of domains. 
What Eq. 5 simply implies is that within the rounding 
width An = g(l*) of n CT (p), the typical domains of the 
vortex liquid of size I* will proliferate in the solid. This is 
because here the system is nearer to the transition point 
of the clean system. For the dimension of interest in the 
present paper, D = 2, we give in Fig. 1 a schematic sketch 
of 5 (0. 

The basic point of observation, which is also one of the 
central points in the present paper, is that g{l) rises from 
zero at I = lo and goes back to zero for I — > oo, and so is 
peaked at the point of inflection 1 = 1*. Consequently, as 
we increase the density of the vortices, n, by increasing 
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FIG. 1. A schematic sketch of the function g(l) for D = 2. 
I = I* is a point of inflection, which is a maximum. 

H, domains of vortex liquid of size I* can first appear 
sharply (sharpness depends on the width of the peak of 
g(l) at I*, which in turn depends on the parameters of 
the system such as a and f(p)) at n — n sm = n cr (p) — 
An = n CT {p) — g(l*)- Also, if £ cr , the finite correlation 
length of the would-be first order transition, is smaller 
than I*, then the domain-size at the sharp crossover is not 
bounded by £ cr , rather it is I*, which can be much bigger 
than £ cr again depending on the system parameters. As n 
is increased more, the function g(l) also satisfies Eq. 3 for 
values of I other than P, and, in particular, the domain 
sizes can be arbitrarily large as An — * 0. The typical 
domain size, however, will always be /*, since by creating 
domains of that size the system gains maximum energy. 

A domain of the vortex solid phase is almost a super- 
conductor, since the vortices are immobile. Similarly, a 
domain of the vortex liquid phase is almost an insulator. 
Low-temperature transport properties of a phase with 
inhomogenous admixture of such domains have been re- 
cently studied in [3]. Because of quantum tunneling of 
the vortices among the domains of vortex liquid, finite 
resistance arises in the system. The formula given in [3] 
for this residual saturation tunneling resistance produces 
good fits to the experimental data [2,4]. One of our aims 
in this paper has been to clarify the origin of such an 
inhomogenous state and to show that it can first occur 
sharply with increasing H. For values of H less than the 
crossover value, H sull the system is a strong supercon- 
ductor with no appreciable residual resistance. This is 
because the domains of the vortex solid are large (the 
domains of the vortex liquid are rare), which makes the 
tunneling of the vortices impossible. For H > H sm the 
system has domains of typical size I* and the resulting 
weak superconductor has low-temperature tunneling re- 
sistance. Because of the proliferation of the domains, 
hysteresis may be observed in the vicinity of H sm , and 



the sharp crossover may look like a first order supercon- 
ductor to metal quantum phase transition. This is the 
remnant effect alluded to before. 

Experiments of [1,2,4] have observed approximate scal- 
ing of resistivity inside the metallic-like phase in the 
phase diagram. This is indicative of a continuous phase 
transition at some critical value of H = H c > H sm . This 
transition, however, must be defined in the sense of per- 
colation, since the vortices are not in the perfect solid 
phase. As we have seen, for H > H sm the sample is in 
an inhomogenous admixture of domains of almost super- 
conductors and almost insulators of typical size I*. As 
H increases towards H a the insulating domains prolifer- 
ate since the energy gain by creating them increases. The 
requisite continuous phase transition, then, can be that of 
two-dimensional randomly diluted Josephson Junctions 
[24] . The fact that the I — V curves are nonlinear at rel- 
atively low currents attests to the presence of Josephson 
Junctions in this region of the phase-diagram [6]. Ig- 
noring capacitive charging energies of the individual do- 
mains, which is reasonable for large enough domain-sizes, 
this system is isomorphic to randomly diluted XY model 
in two dimensions [25]. 

These models exhibit classical phase transitions. When 
the sites or the bonds are randomly diluted, the transi- 
tion temperature T c (x) decreases monotonically with the 
concentration x of occupied sites (or bonds). At the per- 
colation threshold x = x c , which corresponds to H = H c 
in the present problem, an infinite connected cluster of 
occupied sites ceases to exist, and T c (x) = 0. The point 
P : x = x c (H = H c ), T = is a multicritical point where 
long range thermodynamic and geometric order develop 
simultaneously [26]. This true critical point H c may not 
be the same as the would-be first order melting point of 
the clean system, H cr . It will depend on the value of the 
percolation threshold in a more microscopic calculation. 
The experimental observation that H c is about 20 times 
bigger than H sm is a quantitative issue and is beyond the 
scope of the present paper. A numerical simulation study 
of the vortex system to know the function g(l), which in- 
volves parameters C (~ 2ir for a circle), a and f(p) is left 
for future study. Once maximized, this function will give 
an estimate of H sm (n sm = n cr (p) — g(l*)), and will also 
test whether g(l) is indeed a sharply peaked function to 
validate the conjecture of a sharp crossover. In Fig. 2, 
we show a proposed phase diagram of the system within 
the randomly diluted XY model. 

In the conventional theory of the SIT [8] , the Ohmic 
resistivity, p, was predicted to attain a universal value at 
the transition point H c . Its finite temperature scaling- 
form was predicted to be 



where pq is the quantum of resistance, v and z are 
the correlation length exponent and the dynamic expo- 
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FIG. 2. A schematic phase diagram of the system within 
the diluted XY model. H c is the true critical point, close to 
which scaling should apply. H cr is the avoided first order melt- 
ing point. The sharp crossover happens at H sln , beyond which 
domains of size I* proliferate. T C (H) = T* is the transition 
temperature of the diluted XV-niodel, but only a crossover 
temperature for the real two-dimensional system. Below this 
temperature, in a so-called renormalized-classical phase, the 
properties are dominated by dilute thermal excitations above 
the superconducting ground state. 

ncnt, respectively, and F is a scaling function. A small 
non-zero temperature reduces the effective dimension of 
the system from three to two, and a finite-size scaling 
argument produces the scaling-form close to the zero- 
temperature critical point [27] . Recent experiments have 
not found a universal critical resistivity, although the 
scaling form has been verified for T > lOOmK with a 
value of zv ~ 4/3 [1,4,6]. At lower temperatures, break- 
down of the scaling has been observed [4]. 

In our model of the SIT, because P is a continuous 
critical point where correlation length diverges, scaling 
of observables at T — is immediate and follows from 
the well-studied theory of percolation [28]. Scaling at 
finite temperature close to the percolation threshold is, 
however, not well-understood. This is because the zero- 
temperature transition is essentially geometric, and the 
concepts of reduced dimensionality with temperature and 
associated finite-size scaling don't apply. Here, we will 
speculate on an explaination of the experimental observa- 
tions. For any T > at x = x c , p is nonzero since the su- 
perconducting correlation length is finite and is expected 
to be thermally activated. Close to the point P there are 
two relevant correlation lengths that diverge; the perco- 
lation correlation length £ x diverges as £ x ~ \x — x c \~ Ux , 
and the thermal correlation length of the diluted XY 
model, £r, diverges as £t ~ T~ VT [28]. This is true as 
long as £ x > £r, because of the reduced dimensionality 
of the percolating cluster. As this is a crossover effect for 
any x — > x c , but not exactly at the percolation threshold, 



the condition on £p puts a lower bound on temperature, 
only above which the following holds. Here v x and vt are 
critical exponents with values in two dimensions 4/3 and 
0.98— 1.03 [28], respectively. In a scaled expression for p 
for T — > and x —> x c , these two lengths must occur in 
the dimensionless ratio S,t/£,x, 

p(\- - *c|,T) ~ £5P/Kt/€.) - ^(j^r)- (7) 

Here p is a dimensionless constant, is the quantum 
of resistance present for dimensional reasons, and / and 
(f) are scaling functions. The argument of <j> is a power of 
the argument of /. Writing H for x and using the values 
for vt and v x appropriate for 2D, we get 

p(\H-H c \,T)~ ^i}JL^EA). (8) 

in the scaling region close to P. 

In summary, we have proposed a low-field crossover 
from a strong superconductor to an inhomogenous state 
which is a weak superconductor and so can be resis- 
tive even at very low temperature. If the crossover is 
sharp, which is possible in two dimensions but depends 
on the system-parameters, it will look like a first order 
superconductor-metal phase transition. The true ther- 
modynamic superconductor-insulator critical point lies 
within the imhomogenous state and is conjectured to be 
isomorphic to that of diluted XY-model. An explaina- 
tion of the observed finite temperature scaling of resis- 
tivity close to the critical point is presented and that 
it is only approximate, consistent with experiments, is 
pointed out. 
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